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We report here results on the study of the totaUy asymmetric simple exclusion 
processes (TASEP), defined on an open network, consisting of head and tail simple 
chain segments with a double-chain section inserted in-between. Results of numer- 
ical simulations for relatively short chains reveal an interesting new feature of the 
network. When the current through the system takes its maximum value, a simple 
translation of the double-chain section forward or backward along the network, leads 
to a sharp change in the shape of the density profiles in the parallel chains, thus af- 
fecting the total number of cars in that part of the network. In the symmetric case 
of equal injection and ejection rates a = /? > 1/2 and equal lengths of the head and 
tail sections, the density profiles in the two parallel chains are almost linear, char- 
acteristic for the coexistence line (shock phase). Upon moving the section forward 
(backward), their shape changes to the one typical for the high (low) density phases 
of a simple chain. The total bulk density of cars in a section with a large number of 
parallel chains is evaluated too. The observed eff^ect might have interesting implica- 
tions for the traffic fiow control as well as for biological transport processes in living 
cells. An explanation of this phenomenon is offered in terms of finite-size dependence 
of the effective injection and ejection rates at the ends of the double-chain section. 

Keywords: TASEP; traffic flow models; non-equilibrium phase transitions; traffic 
on networks; biological transport processes 
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I. INTRODUCTION 

Traffic often takes place along linear tracks which are interconnected to form a network 
structure. However, in the presence of hard-core exclusion the collective behavior of such 
transportation networks is not yet well understood. Indeed, the role of junctions has been 
considered only recently on the example of simple networks involving no more than two 
junctions, see [H [2] or on self-similar tree-like topologies, see |3j and references therein. 

The idea of studying networks composed of chain segments, which exhibit the bulk be- 
havior of an open TASEP under boundary conditions, given in terms of effective input and 
output rates, was first advanced in our work [T]. The network considered there consists of 
two vertices of degree 3: one of out-degree 2 and the other one of out-degree 1, connected by 

2 chains of the same direction; the third edge of each of these vertices belongs to a directed 
chain coupled to a separate reservoir of particles, see Fig. [Tj The appearance of correlation 
effects close to the ends of the chain segments, as well as of cross-correlations in the double- 
chain segment was demonstrated. The same approach was applied in Ref. P] to an open 
network consisting of one vertex of degree 3 and out-degree 1. The two incoming chains 
are coupled to one reservoir, and the outgoing one is coupled to another reservoir. Different 
versions of simple networks were studied in Refs. |1] and |5]. In the former reference two 
cases were investigated: (a) two vertices of degree 3 connected by three chains, one of which 
has the opposite direction to the remaining two (closed system); (b) two vertices of degree 

3 connected by two chains with the same direction; the remaining incoming and outgoing 
chains are coupled to reservoirs with the same particle density (open counterpart). In Ref. 
[5] graphs containing vertices of degree 4 and out-degree 1, 2, and 3 were considered. The 
notion of particle-hole symmetry in the presence of a junction was carefully analyzed and 
an appropriate interpretation on the microscopic level was given. TASEP with parallel up- 
date on single multiple-input-single-output junctions has been investigated too ^ . Clearly, 
the above works have treated TASEP on diverse, but simple fixed network topologies. The 
main concern was the construction of the phase diagram under different open boundary 
conditions. 

Complex networks have also been a focus of research in the last decade. Attention 
has been paid to the traffic fluctuation problem in networks: the dependence between the 
mean value of the of traffic passing through a node (or a link) in a time interval and its 
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FIG. 1: (Color online) Schematic presentation of the considered system: TASEP, defined on an 
open network (solid black line), consisting of head Ci and tail C4 simple chain segments with a 
double-chain section (C2 and C3) inserted in-between. In the system, considered in [1] all the 
simple chain segments have equal lengths Li = L2 = L3 = L4 (more details are in the text). The 
light grey (green) dash-dotted line shows a version of the network, where the double-chain section 
has been translated forward by AL so that the head and tail simple chain segments have lengths: 
Li = Li + AL and L4 = L4 — AL. 

standard deviation has been studied. It has been recognized that in the evolving networks 
functional units may emerge, which may be represented by topologically distinct subgraphs, 
see [7] and references therein. Different types of diffusive dynamics, like spreading of disease 
and traffic or navigated walks have been studied on different networks, see, e.g., p-[TU]. 
In particular, high-density traffic of information packets on sparse modular networks with 
scale-free subgraphs was studied, see [H]. Most of the results of graph theory relevant to 
large complex networks were related to the simplest models of random graphs. A large 
amount of the research was devoted to scale- free networks, i.e., networks with power-law 
vertex degree distribution. Traffic rules on such graphs usually include particle creation 
at randomly selected nodes and mutually interacting random walks to different specified 
destinations, see |12]. In [13] the authors studied the stationary transport properties of 
the TASEP with random-sequential update on complex networks, both deterministic and 
stochastic. The TASEP rules were generalized, by an obvious extension of the rules applied 
in [1] , to fixed connected networks of Ng directed segments (each consisting of L ^ 1 sites) 
and nodes. The theoretical approach was based on a combination of models of complex 
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networks with the well-known mean field (MF) results for the simple chains ('segments') 
which are assumed to connect the vertices of the underlying directed graph. These chain 
segments, representing the edges of the network graph, have to be long enough to make 
reasonable the application of the MF results. At that, the correlations that may build up 
close to the nodes of the network have to be neglected. 

Recently, applications to biological transport have motivated generalizations of the 
TASEP to cases when the entry rate is chosen to depend on the number of particles in 
the reservoir (TASEP with finite resources) [13]. Last year, the cases of multiple competing 
TASEPs with a shared reservoir of particles [15] and with limited reservoirs of particles and 
fuel carriers [H] were studied too. 

II. MODEL AND NUMERICAL SIMULATION RESULTS 

In the asymmetric simple exclusion process, hard-core particles move along one- 
dimensional lattice of sites, which can be occupied by one particle at most {tj = 0, 1}, i = 
1, . . . ,L. The particle can move right (left) with probability p (1 — p) to a nearest neighbor 
site, if the site is empty. In the extremely asymmetric case particles are allowed to move 
in one direction only — this is the totally asymmetric simple exclusion process (TASEP). 
Its steady states are exactly known for both open and periodic boundary conditions, for 
continuous-time and several kinds of discrete-time dynamics. Here, we focus our attention 
on the steady states of the open TASEP with continuous-time stochastic dynamics, mod- 
eled by the so called random-sequential update. For a review on the exact results for the 
stationary states of TASEP, under different kinds of stochastic dynamics, and its numerous 
applications, we refer the reader to [T71 - [T^ . 

Our goal here is to present some interesting new effects, observed in a TASEP, defined 
on a simple network, consisting of head and tail simple chain segments (Ci and C4) with 
a double-chain section (C2 and C3) inserted in-between. The model system is presented 
schematically in Fig. [T| The particle injection rate at the left end of the network is a and 
the particle ejection rate at the right end of the network is /3. Pi, is the branching point 
— the last site {i = Li) of the head chain segment Ci, where particles can take with equal 
probability pj the upper C2 or the lower C3 branch of the double-chain section. Pm is the 
merging point (first site of the tail segment i = Li + L2 + 1), where the particles moving 
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FIG. 2: (Color online) Density distribution p(x) (from numerical simulations) along the different 
chain segments of the open network of the type, shown in Fig. 1, when a = 0.75, /? = 0.75. The 
solid (red) squares are the reference results for the network with the (MC,CL,MC) phase structure 
appearing when Li = L2 = L3 = L4. Simple translation of the double-chain section, while the 
lengths of the double chain segments C2 and C3 are kept fixed L2 = L3 = 50, and also the length of 
the whole network is kept fixed at L = 150 sites, causes a noticeable change of the density profiles 
of C2 and C3. Translation forward leads to a density distribution shape on C2 and C3 (shown 
with solid (green) down-triangles), characteristic of the LD phase, while translation backward to 
a density distribution (solid (blue) up-triangles), characteristic of the HD phase. For comparison 
the density distribution (in a MC phase) of a simple chain of length 150 sites is also shown with 
empty grey circles of larger size. 

along the C2 and C3 chain segments merge. The particle current in the system is denoted 
by J (it is in the same direction for both chains of the double-chain section). 

The phase structure of the network, studied in [I] when all the simple chain segments 
have equal length, L\ = L2 = Lj, = L4, is presented by the triplet (Xi, ^2,3, X4), where 
Xn in = 1,2,3,4) stands for one of the stationary phases of the simple chain segment Cn'- 
LD — low density, HD — high density, MC — maximum current, and CL — coexistence 
line. Our analytical analysis of the allowed phase structures, based on the properties of 
single chains in the thermodynamic limit, and the neglect of the pair correlations between 
the nearest-neighbor occupation numbers at the junctions of the different chain segments, 
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FIG. 3: (Color online) Nearest-neighbor correlations Fcorr in the network with a double-chain 
section as a function of the scaled distance x for a = 0.75, /? = 0.75. The simulation results for 
Fcorr correspond to the density distributions, shown in Fig. 2, for the different positions of the 
double-chain section. 

yielded 8 possibilities. Here we focus our investigation on 3 of the most interesting cases 
(MC,LD,MC), (MC,CL,MC), and (MC,HD,MC) which appear when the boundary rates 
satisfy the inequalities a > 1/2, (3 > 1/2, corresponding to the maximum current phase of 
a single chain. We have shown that the phase state of the chains in the double-chain section 
depends on the effective injection rate a* of particles at the first site of each of the chain 
segments 6*2,3 and on the effective removal rate /3* of particles from the last site of each of 
these chains. Our further studies of the system reveal a rather interesting property whenever 
the network with Li = L2 = = L4 has the (MC,CL,MC) phase structure. One can change 
effectively the effective rates a* and /3* at the network junction points simply by changing 
the position of the double chain section along the network, while keeping fixed both: the 
total length L of the network and the length L2 = L3 = Ld of the double-chain section. This 
is clearly observable in Figs. [2]-|4j As one can see in Fig. [2| a simple translation forward 
of the double chain section leads to a density profile on C2 and C3 characteristic of the LD 
phase (shown with solid down-triangles), while translation backward induces a density profile 
characteristic of the HD phase (shown with solid up-triangles) . As a reference, the results 
for the density profiles of the system with segments of equal length Li = ^2,3 = L4 = 50 are 
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FIG. 4: (Color online) Density distribution p{x) (simulation results) along the two equivalent 
chain segments of the double chain section for bigger size network (total length L = 300 sites) 
at a = 0.75, /3 = 0.75. The results show the increased sensitivity of the density distributions to 
even a very small position change AL of the double chain section. The solid (red) squares are the 
reference results for Li = L2 = L3 = = 100 (AL=0). 

shown with solid squares. 

The characteristic change of the density profiles is accompanied also with the correspond- 
ing characteristic change in the nearest-neighbor correlation function F^orr, displayed in Fig. 
|3} for the same cases, as the ones shown in Fig. [2j For comparison the density distribution 
(in the MC phase) and nearest-neighbor correlations of a simple chain of length L = 150 
sites are also shown in both figures with empty grey circles of larger size. The Monte Carlo 
simulation results, presented here, are for a relatively small system of fixed total length 
L = Li + ^2,3 + L4 = 150 sites and fixed size of the double-chain section, L2 = = 50. The 
ensemble averaging was performed over 200 independent runs and 1.5 x 10^ Monte Carlo 
steps were omitted in order to ensure that the system has reached a stationary state. 

Numerical study of a larger system with total length L = 300 sites reveals even higher 
sensitivity of the density distribution p{x) along the chains C2 and C3 of the double-chain 
section with respect to quite small changes of the loop position on the network, see Fig. |4} 
Already at AL = 2 one can observe a noticeable change in the density distribution. 
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III. THEORETICAL ANALYSIS 



A. Mean field theory 



The possible phase structure of the present network was analyzed in our paper by 
using the exact results for simple bulk chains C„, n = 1,2,3,4, coupled to each other by 



means of effective injection and removal rates. To be specific, let t- 



in) 



1 9 T 

±, z,, . . . , ±jni 



denote the particle occupation number of site i in chain C„ which contains L„ sites. Here we 
let the length of the double-chain defect L2 = L-^ = and the total length of the network 
Li + + L4 = L to be fixed. The mean field theory assumes that all L„, n = 1,2,3,4, 
are large enough, so that the average occupation number {t^'^^){L; a, P) of site i G C„ in 
the stationary state of the whole system, under injection rate a and removal rate (3, can be 
identified with the average occupation number {T^'^^){Ln', a*, /?*) in the stationary state of a 
single chain of length L„ under some effective injection (a*) and removal (/?*) rates. Note 
that at the open boundaries of the whole network the continuity of the current Jl through 
the system implies 



{4)\ 



(1) 



Under the above assumptions, these equalities will be approximated by 



(2) 



Similarly, by neglecting the nearest-neighbor correlations at the junctions of different chains, 
the current continuity equation yields the approximate equalities 



Jl ^ (ri^^))(Li;a,/3*) \ - 1- {rr){L,; a*,, 



_(2)\ 



-(3) 



.(4)\ 



(3) 



Hence, taking into account the equivalence of chains C2 and C3, we set 0.2 

f^l = 1^1 = l^*d-: define 



a 



di 



/3r = l-(rP))(L,;aS,/3*), 



1 



rf))(L4;a:,/3) 



(4) 
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or, alternatively, 

a J = jTT-TTs , a A = 7TT . (o) 

2[l-(rf'^))(L,;a*,/3*)] 1 - {ri'^){L,; a*,, P) 

Thus, definitions ^ express the effective injection (removal) rates of chains €2,3, C4 (Ci, 
C2,3) in terms of the average occupation number of the last (first) site of the preceding (sub- 
sequent) chain, while definitions ([S]) express the corresponding effective injection (removal) 
rates in terms of the finite-size current and the average occupation number of the first (last) 
site of the same chain. The consistency of these definitions is a measure of the extent to 
which nearest-neighbor correlations at the junctions can be neglected. 

The above equations essentially simplify when L„ ^ 1, n = 1,2,3,4, by replacing the 
finite-size properties of the current and the density profiles with their thermodynamic coun- 
terparts. To this end we define 

Pin\(^^P) = lim (r}"^)(L„;a,/3), pi"J(a,/3) = lim (rj^"^)(L„; a, /3), 
pi:L(«,/3) = lim -5^(r("))(L„;a,/3). (6) 

L,i-s>oo Ln — 
i=l 

In this case the current is determined by the bulk densities (omitting the dependence on the 
injection and removal rates): 

J = lim Jl = piifi(l - piifi) = 2pS(l - pS). (7) 

Hence, we find expressions for the possible bulk densities in terms of the current in the 
thermodynamic limit: 

pZ = P±iJ), pil = P±iJ), pS=P±(-//2), (8) 

where 

p±(j) = ^(i±v/r^). (9) 

In the most interesting case a > 1/2 and /3 > 1/2, both chains Ci and C4 are in the 
maximum current (MC) phase [Ij, when the exact results in the limit L — 00 yield 

plJ(a:,/3) = l--^, pJ,:^,(a:,/3) = -^. (10) 



(MC) 


= pil(MC) = 


1/2, 


1 

' 4^' 


piiU«,/3i*) = 


1 

4^' 


1 

4a:' 


pS(«:,/3) = 


1 
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In this case from Eq. (8) it follows that either p\^i^. = p_(l/8) ~ 0.146 and the defect chains 
C2,3 are in the low density (LD) phase, or p^^^f^ = p+(l/8) ~ 0.854 and the defect chains 
6*2,3 are in the high density (HD) phase. 
Consider first the LD case, when 

pSf /3d) = pS(«d, p:) = p£f («d, p:) = «d(i - «d)/3d- (11) 

Then, in the large L limit Eqs. Q and ^ yield 

Consider next the HD case, when 

P?n'\<, /3d) = 1 - /3d(l - /3d)/«d, P£f /3d) = pS(«d, /?:) = 1 - /3d- (13) 
Then, the large L limit of Eqs. Q and ([s]) yields 

«:=^pi'i(«,/3i) = ^, a: = ^=2p+(l/8), /5^ = l-pS(a:,/3:)=p_(l/8). (14) 
We see that in both cases there remain undetermined effective rates: in the low density 



case these are j3} and the related al, see Eq. (12), and in the high density case and 



the related /3l, see Eq. (14). We explain this feature by the fact that the relation between 
the magnitudes of and (3^ depends on the phase of the chains C2,3 when chains Ci 
and C4 are in the maximum current phase. Then the current through the whole system 
attains its maximum value J = 1/4 and the bulk densities in Ci and C4 take the only 
possible value pli^Jn^ = p[^ii, = 1/2. However, the value of the current through each of the 
chains 6*2,3, •^''2,3) = J/2 = 1/8 allows them to be either in the low density phase, when 
Pd ^ — P-(l/8)5 or in the high density phase, when > f3^ = p_(l/8), or on the 
coexistence line, when = = p_(l/8). This means that there is a degree of freedom 
due to the fact that the average density of particles in the chains C2,3 is not fixed. In the 
next section we show how one can deduce the missing injection/ejection rate from the fit of 
the density profile with the exponential distribution predicted by the domain wall theory. 

Our former computer simulations [IJ have shown that in the symmetric case of Li = 
L4 = 200 and a = P > 1/2 the double-chain segment is found to be on the coexistence 
line = < 1/2 (known also as a shock phase) when a completely delocalized domain 
wall exists. In this case the local density profile is linear and changes in the interval from 
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Pb„ifc(LD) to pbuik{}iD), so that prf(a^,a^) = 1/2. Evidently, for a finite-size system the 
effective rates depend on the lengths of all the chains, Li, and L4. Simple arguments 
lead to the conclusion that as Li — )■ 1, the effective injection rate monotonically increases, 
so that —> a/2 > 1/4, while (32 = p-(l/8) < 1/4 does not change significantly if ^ 1, 
the defect chains 6*2,3 go into the high density phase. On the other hand, when L4 — )■ 1, 
monotonically increases, so that — )■ /3 > 1/2, while = p_(l/8) < 1/4 does not change 
significantly, the defect chains 6*2,3 go into the low density phase. That is actually observed 
in our present computer simulations, the results of which are illustrated in Figs. |2} The 
noticeable deviations of the density profiles from the standard infinite chain counterparts 
are due to finite-size effects, as well as to the correlations that appear near the points of 
inhomogeneity of the network, see Fig. [3j The sensitivity of the effect increases with the 
length of the chains, see Fig. |4], since then the smeared phase transition in 6*2,3 becomes 
sharper. 

Some numerical data for a network with fixed L = 150, = 50 and changing Li = 
100 — L4, are given in Table 1. The values of and are calculated from the finite-size 
numerical data for the local densities (t^^), (t^^''), {t[^'^^) and (r^f^), by using equations 
Q and Q with Jl = J(MC) = 1/4. 













from Eq. (4) 


from Eq. ( 5 ) 


from Eq. (4) 


/3d 

from Eq. (5) 


25 


0.472 


0.856 


0.436 


0.856 


0.236 


0.222 


0.144 


0.146 


50 


0.299 


0.824 


0.175 


0.791 


0.150 


0.152 


0.176 


0.158 


75 


0.281 


0.590 


0.144 


0.429 


0.140 


0.146 


0.410 


0.291 



Table 1. 



As is seen, the relative values of the effective rates change with the defect position, 
according to our expectations, and properly describe the phase changes, observed in the 
density profiles. The largest numerical discrepancies are observed between the values of 
obtained from Eqs. Q and ([s]) in the cases of networks with Li = 50 and Li = 75. Probably, 
they can be explained by taking into account the high nearest-neighbor correlations between 



12 



the defect section and the tail chain. 



B. Domain wall theory 

According to the Domain Wall (DW) theory, the configurations of TASEP on a simple 
chain with open boundaries can be approximated by two regions, of low and high density, 
separated by a domain law of zero width [20|, see also [21j. Thus, all the configurations of a 
chain of L sites can be labeled by a single integer fc = 0,l,...,L, such that sites < i < /c 
belong to the low density phase with uniform density p_, and sites k + 1 < i < L belong to 
the high density phase with uniform density p+. The extremal values of A; = and k = L are 
understood to label configurations corresponding to the pure high and low density phases. 

We recall that the probability P{k,t) for finding at time t the domain wall at position k 
satisfies a master equation with refiecting boundary conditions at = and k = L. The 
corresponding stationary solution P*(/c) has the simple exponential form 

P,(A;) =r-VAr, (15) 

where 

a{l-a) |exp(-l/0, « < /? < 1/2 
- /5) [ exp(l/0, /3 < a < 1/2 ' 
where ^ > is the domain wall localization length; the normalization factor A/" is 

^f■■=J:--' = ^^^. (16) 

k=0 

Thus, the local density of particles p{i) := (rii) is given by 

p(0 = p^N-'Y^r-' + p_A/-i X^r-'^ = p„ + (p+ - p-) ^^_~!2i - (17) 

fc=0 k=i 

This expression describes different shapes depending on the value of r. 



(a) On the coexistence line (CL) r = 1 and Eq. (17) yields the linear profile 



pcLii)=P- + {p+-p-)j^, l<i<L. (18) 

(b) In the low density (LD) phase r < 1 and up to exponentially small corrections one 
obtains 

Pld{t)= P- + {p+- P-)r''-'^\ r<l, l<i<L. (19) 
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(c) In the high density phase r > 1 and up to exponentially small corrections one obtains 

PHoii) =P+- {p+ - P-)r-\ r > 1, l<i<L. (20) 

We shall concentrate on the last two cases which will be considered on the macroscopic 
scale i/L = X. To this end it is convenient to introduce the parameter t = ^/L and rewrite 



expression (19) as 

Pld(x) =p_ +Aexp[(x- 1- l/L)/t], L-^<x<l, (21) 



where, Pld{.x) = pld{xL) and A = p+ — p_. Similarly, expression (20) takes the form 



p^^(x) = p+ - Aexp(-x/t), L'^<x<l, (22) 

where, Phd{x) = phd{xL) and A = p+ — p_. 

The results from the interpretation of our numerical results on the local density profiles 
within the DW theory are explained in the following two subsections. 

1. Interpretation of the density profile fit in the LD phase 

The fit to the simulated local density profile by function (21) is shown in Fig. |5} Its 
very high quality yields reliable estimates of p_ = ~ 0.148, A = 0.421 ± 0.002 and 
t = 0.0488 ± 0.0003, hence, r = exp(-l/0 ^ 0.664. 

The value of = 0.148 is in fairly good agreement with the average current J^'^'^^ = 
jmax 12 through each of the chains 6*2,3: ajj(l — a*^ ^ 0.126; the difference between 0.126 
and the expected value j'^"^/2 = 1/8 through each of the chains of the double-chain section 
can be attributed to finite-size effects in a network of total length L = 150. 

The most important observation is that, by using the prediction of the DW theory, 

(3*,{l-P:)=a:{l-a:)/r = j(''^^/r, (23) 

we can evaluate the effective ejection rate (3^ by solving the above equation at r = 0.664, 
with the result 

/3* = (1/2) [1 - -^1^47(2^] - 0.255. (24) 

Note that this value satisfies the conditions for the existence of a LD phase, but differs 
significantly from the mean field estimates given in Table 1. To check the validity of the DW 
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FIG. 5: (Color online) The simulation results for the particle density profile p{x), x = {i — Li)/ L2, 
i = Li + 1, Li + 2, . . . , Li + L2, in the LD phase of a chain C2,3 in the network, are shown by 
centered grey circles, the fit with the function ( |21[ ) is shown by solid (red) line. The quality of the 
fit is characterized by the statistical criteria /DoF = 1.0004 x 10~^ and = 0.99963. 

theory, we have performed computer simulations of TASEP on a single chain with injection 
and ejection rates a = 0.148 and (3 = 0.255, respectively. The corresponding density profiles 
are compared in Fig. |6| One sees a qualitatively very close behavior of the density profiles, 
although in the region of upward bending, close to the right-hand end of the chains, the local 
density at the single chain sites is considerably higher than the one at the corresponding sites 
of the C2,3 branches. At the right end this difference reaches its maximum: p^'^(L2) — 0.429, 
while for the simple chain p^^{L2) — 0.495. The latter value is not influenced significantly 
by finite-size effects, because it coincides, within numerical accuracy, with the value of the 
density in the thermodynamic limit for a simple chain with the same boundary rates: 



We explain the above discrepancy by the presence of nearest neighbor correlations between 
the last site of the chains €2,2, and the first site of the tail chain C4: 



From the exact expression for the current J^^'^^ through each of the chains 6*2,3 if follows 




(25) 




(26) 
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FIG. 6: (Color online) Comparison of the simulation results for the particle density profile p{i), 
i = 1, 2, . . . , 50, in the LD phase of a chain C2,3 in the network (shown by solid (red) circles) and 
an open simple chain (shown by solid (blue) squares) under the estimated rates a = 0.148 and 
P = 0.255. 



that 

7(2,3) I I?('^'4) 

= ' ■ (27) 

1 -Pl 

Hence, by inserting the numerically evaluated Fm^^ — 0.048 and p'f' ^ 0.590, we obtain 
pf^'^^ ~ 0.424 which is fairly close to the simulations result p'^'^{Ld) ^ 0.429. 

Finally we note that although the DW theory gives an excellent description of the density 
profile of chains in the LD phase, the maximum reached at the right-hand end of a long 
simple chain is lower than the bulk density p+ = 1 — of the high-density phase that 
supports the same current for all a*^ < f3^: 

hm pUL) = - < (1 - «d) = P+- (28) 

Indeed, the interpretation A = — p^, given by the DW theory, yields a value for the 
bulk density of the HD phase, p+ = A + p_ ~ 0.569, which cannot support the current 

j(2,3) jmax jc^ 
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FIG. 7: (Color online) The simulation results for the particle density profile p{x), x = {i — Li)/ L2, 
i = Li + 1, Li + 2, . . . , Li + L2, in the HD phase of a chain C2 3 in the network, are shown by 



centered grey circles, the fit with the function (22) is shown by solid (red) solid line. The quality 



of the fit is characterized by the statistical criteria / DoF = 9.134 x 10 and E? = 0.99986. 



2. Interpretation of the fit in the HD phase 



The fit to the simulated local density profile by function (22) is shown in Fig. [7| Its 
excellent quality yields rehable estimates of p+ = 1 - ~ 0.853 > 1/2, A = 0.594 ± 0.002 
and t = 0.0570 ± 0.0002, hence, r = exp(l/Lt) ~ 1.42. 

From the relation for the bulk density in the HD phase, p+ = 1 — we deduce the value 

^ 0.147 for the effective ejection rate. This gives a very good estimate of the current 
j(2,3) ^ _ ^ 0.1254 through each of the chains ^2,3- 

Next, following the DW theory, we solve the equation 



a, 



at r = 1.42, to obtain 

al = (1/2) [1 - Vl - 4r J(2.3)] ~ 0.232. 
This value agrees fairly well with the mean field estimates given in Table 1. 



(29) 



(30) 



To check the consistency of the fit, from Eq. (20) we obtain 



p_ ~ Phd{1) =P+- A/r ~ 0.435, 



(31) 
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FIG. 8: (Color online) Comparison of the simulation results for the particle density profile p{i), 
i = 1,2, . . . , 50, in the HD phase of a chain C2,3 in the network (shown by solid (red) circles) and 
an open simple chain (shown by solid (blue) squares) under the estimated rates a = 0.232 and 
P = 0.147. 

which almost coincides with the numerical value of 0.436 for the density at the first site of 
the chains 6*2,3. On the other hand, the analytical expression for the density profile in the 
thermodynamic limit leads to the estimate 

lim phd{1) = 1 - /3^(1 - P:)/a: - 0.459, (32) 

which coincides, within numerical accuracy, with the numerical value PHi){l) — 0.460 ob- 
tained for a simple chain under the same boundary rates, see Fig. |8} and which is about 5 
percent higher than the numerically obtained value phd{^) = 0.436. This small discrepancy 
is due to the small nearest-neighbor correlations between the last site of the head chain and 
the first site of chains 6*2,3 : 

FiY^ = (r«rf '^^) - (r«) (rf •^)) =^ 0.011. (33) 

To check the validity of the DW predictions, we have simulated the density profiles in a 
single chain with L = 50, a = 0.232 and /3 = 0.147. A comparison of the results is given in 
Fig. i 

However, as in the previous case, the DW interpretation A = p^~p^ leads to p_ ^ 0.259, 
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a value which does not support the current J*^^'^^ = J^°-^ and, therefore cannot describe a 
bulk LD phase. 



Non-equilibrium phenomena are much more often encountered and more diverse in nature 
than the true equilibrium phenomena. The development of our understanding of physics far 
from equilibrium is currently under way. In this respect the study of simple non-equilibrium 
models plays very important role. The asymmetric simple exclusion process is one of the 
simplest non-equilibrium models of many-particle systems with particle conserving stochastic 
dynamics and boundary induced phase transitions. 

Here, we have reported results on a rather unexpected property of the TASEP, defined 
on a simple network consisting of a single chain with a double-chain insertion. In the 
symmetric case of equal lengths of all the segments, and at equal injection and ejection 
rates a = (3 > 1/2, the network has the maximum-current (MC,CL,MC) phase structure. 
Then, a simple translation of the double-chain section forward or backward along the single 
chain induces a qualitative change in the local density profiles in the parallel chains which 
is characteristic of a phase transition. This "position-induced phase change" is caused by 
the change of the effective rates a* and /3* at the two junction points of the network. 

Using the continuity of the finite-size current Jl one can determine the effective rates 
a* and /3* at the two ends of the loop. Our theoretical analysis, based on the results for 
infinitely long chains, leads to values of a* and /3* which change with the double-section 
position and describe well the phase changes obtained numerically in the profiles. Note, for 
a finite-size system the effective rates depend on the lengths of all the chains, Li, Ld and 

u. 

Quite interestingly, upon increasing the number of parallel chains in the inserted section, 
provided the parallel chains are in the LD phase, the total bulk density of particles in them 
goes down to to the value 1/4. Indeed, in the case of n parallel equivalent chains the current 
through each of them is 3'^°'^ jn = l/(4n). The low density phase that supports this current 
has a bulk density, compare with Eq. Q, 



IV. DISCUSSION 




(34) 
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Therefore, 

lim np^{J""'^/n) = lim n - f 1 - Vl - n-A = 1/4. (35) 

At that, the total bulk density in the n parallel chains, when they are in the HD phase, 
increases linearly with n 

np+{J"'^yn) = n - 1/4 + 0{n-^), as n ^ oo. (36) 

Another novel observation, which we cannot explain completely, concerns the validity of 
the DW description of the observed density profiles in the LD and HD phases of chains 6*2,3 . 
Apparently, see Figs. |5] and [7| these profiles are excellently fitted by a single exponential 
function, which implies the existence of a well-defined localization length. However, the 
value of the amplitude A of the exponential function does not agree with the DW prediction 
A = pj^ — p^, where p+ and p_ are the bulk densities of the HD and LD phases, respectively, 
which support the required current J'""^/2. This phenomenon cannot be attributed to 
correlations specific to the network under consideration, because the same conclusions hold 
true for the simple chains under the corresponding boundary rates. 
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